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Abstract 

We examine a multidimensional optimization problem in the trop- 
ical mathematics setting. The problem is to minimize a nonlinear 
function defined on a finite-dimensional semimodule over an idem- 
potent semiheld under linear inequality constraints. Given two arbi- 
trary matrices, the objective function is defined by one of the matrices 
through the use of multiplicative conjugate transposition, whereas the 
inequality constraints are given by the other matrix. We start with 
an overview of known tropical optimization problems with both linear 
and nonlinear objective functions, subject to linear equality and in- 
equality constraints. A short introduction to tropical algebra is given 
to provide a formal framework to solve the optimization problem under 
study. As a preliminary result, a general solution to a linear inequal- 
ity with an arbitrary matrix is presented. We describe an example 
optimization problem drawn from project scheduling and then offer a 
general representation of the problem. Furthermore, we introduce an 
additional unknown variable and reduce the problem to solving a lin- 
ear inequality, where the new variable plays the role of a parameter. 
A necessary and sufficient condition for the inequality to hold is used 
to evaluate the parameter, whereas the solution of the inequality is 
taken as a solution for the optimization problem. As a result, under 
fairly general conditions, a complete direct solution to the problem is 
derived in a compact vector form. Numerical and graphical examples 
for two-dimensional problems are given to illustrate obtained solutions. 
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1 Introduction 



Tropical (idempotent) mathematics, as an applied mathematical theory and 
methods based on the notion of idempotent semirings, has its origin in a 
few works, including [D 121 El E]> all inspired by problems from operations 
research. Over the last decades, there have been many theoretical results 
achieved and various applications developed in a number of studies, which 
are reported in such monographs as [SJ El El El El Ell [HI 121 E31 E] , as well 
as in a wide range of contributed papers. 

Even though models and methods of tropical mathematics find expand- 
ing applications in many fields, including engineering, computer science, and 
economics, the scope of operations research remains to be the main domain, 
which supplies more new problems to solve and uses more new results to 
apply. One of the areas of investigation directly concerned with problems in 
operations research is tropical optimization, which deals with analysis and 
solution of optimization problems formulated in the tropical mathematics 
setting. 

A minimax earliness problem in machine scheduling in |15| is an early 
instance of real-life problems that can be represented and solved within the 
framework of tropical mathematics. Further examples include multidimen- 
sional optimization problems appeared in location analysis, transportation 
networks, decision making, and discrete event systems. 

The optimization problems usually take the form of minimizing (max- 
imizing) functions defined on finite-dimensional semimodules over idempo- 
tent semifields, and may have additional equality and inequality constraints 
imposed on the feasible solutions. Some problems involve both an objec- 
tive function and constraints that are linear in the tropical mathematics 
sense. These tropical linear optimization problems, which appear to be for- 
mal analogues of those in conventional linear programming, were apparently 
first considered in [TBI El E] • 

In another important class of problems dating back in the literature to 
the works 17]. the constraints are still linear, whereas the objective func- 
tion is not. The class actually combines problems with functions involving a 
conjugate transposition operator, which does not preserve tropical linearity. 

There are a number of solution techniques developed to handle the tropi- 
cal optimization problems. Some problems have complete solutions obtained 
in an explicit vector form in terms of a general semiring as those in the early 
works [16 1 1151 [TT1 IT]. For other problems, existing solutions are given for the 
case of a particular semifield in the form of iterative algorithms that give 
particular solutions if any, and indicate that there is no solution otherwise 
(see, for instance, [T51 [T5]). 

In this paper, we examine a multidimensional optimization problem that 
appears, in particular, in project management when an optimal schedule is 
constructed to minimize flow time of activities in a project under various 
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activity precedence constraints. 

We formulate the problem in terms of a general semimodule over an 
idempotent linearly ordered radicable semifield. Given two arbitrary matri- 
ces, the problem is to minimize a nonlinear objective function defined by one 
of the matrices through multiplicative conjugate transpose, and involves lin- 
ear inequality constraints given by the other matrix. The problem presents 
a sufficient extension of that considered in [20], where an irreducibility con- 
dition is imposed on the matrices, whereas the inequality constraints are 
less general. Solutions to an unconstrained version of the problem together 
with related applications in location analysis and stochastic discrete event 
systems are presented in [2T1 1221 [1311231121] . 

To solve the problem, we implement and further develop results offered in 
[251 E21 [13] , including solutions to linear inequalities and extremal properties 
of eigenvalues of matrices. We follow the approach proposed in [20], which 
introduces an additional unknown variable and then reduces the problem 
to solving a linear inequality, where the new variable plays the role of a 
parameter. A necessary and sufficient condition for the inequality to hold 
is used to evaluate the parameter, whereas the solution of the inequality is 
taken as a solution of the initial problem. As a result, under fairly general 
conditions, a complete direct solution to the problem is derived in a compact 
vector form, which is suitable to both further analysis and applications. 

The rest of the paper is organized as follows. We start with an overview 
of some known multidimensional tropical optimization problems with both 
linear and nonlinear objective functions and linear constraints. We briefly 
discuss related solution methods and indicate existing areas of application. 

A short introduction to tropical algebra is then given to provide an 
appropriate formal framework in terms of a general idempotent semifield to 
solve the optimization problem under study. The introduction includes clear 
graphical illustrations intended to clarify main definitions and basic facts by 
the example of the classical semifield (R max ,+ • 

Furthermore, a complete explicit solution for a linear inequality is ob- 
tained as an important prerequisite. The solution extends known results for 
the inequality with an irreducible matrix to the case of irreducible matrices. 

Finally, we offer an example problem drawn from project scheduling 
and then give a formal description of the optimization problem of interest. 
We exploit the above solution for inequalities to obtain our main result, 
a complete solution to the problem. Numerical examples of solving two- 
dimensional problems and related graphical illustrations are also included. 

2 Tropical optimization problems 

Multidimensional tropical optimization problems constitute a significant re- 
search domain in tropical mathematics, dating back to the works |16t 1 15 j . 
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The problems arise in real-world applications in various areas. Specifically, 
[El 021 ES1 113 ESI EH] provide topical solutions to problems in project 
scheduling. Furthermore, [301 E31 EH E] solve problems in location analysis. 
Finally, there are applications developed in transportation networks [71 126j. 
in decision making [32], and in discrete event system analysis |21j . 

Below we give a short overview of one class of multidimensional opti- 
mization problems and briefly discuss their solution methods. The prob- 
lems are formulated as to minimize functions defined on vectors of a finite- 
dimensional semimodule over an idempotent linearly ordered radicable semi- 
field. Both unconstrained and constrained problems are considered, where 
the constraints have the form of linear vector equations and inequalities. 

We present optimization problems with both linear and nonlinear ob- 
jective functions. However, only those nonlinear functions are under study 
which are composed with a conjugate transposition operator. As our overview 
shows, a substantial part of problems known in the literature directly falls in 
this category or can be placed into it after some equivalent transformations. 

In the description of the problems, we use the symbols A, B, and C 
for given matrices, b, d, p, and q for given vectors, and x for the unknown 
vector. Matrix and vector operations are considered as defined in terms of an 
appropriate idempotent semimodule. The minus in the exponent denotes the 
vector multiplicative conjugate transpose. Further details on the notation 
used here and throughout the paper, and related definition are given in the 
next section. 

2.1 Linear objective functions 

One of the long-known and well-studied tropical optimization problems is a 
direct tropical analog of linear programming problems, which is to 

... t 

minimize p x, 

subject to Ax > d. 

Complete closed-form solutions to the problem together with related 
duality results are given under various algebraic assumptions in a number of 
publications. Specifically, |16|. IT] examine the problem in the case of general 
semirings; the former offers a solution based on an abstract extension of 
the conventional linear programming duality, and the second suggests a 
residual-based solution technique. Furthermore, [T7] develops results in [TB] 
to concentrate on the idempotent seniifielcL R m ax .+ •> 

whereas [33] offer a 

solution for IR max ,+ and R max , x in the context of the theory of max-separable 
functions. 

A tropical optimization problem in R maXi + with more constraints, which 
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can be written in the form 

... j> 
minimize p x, 

subject to Ax < d, Cx > b, 

g < x < h, 

is considered in |34|, [35j [36[ |26] within the framework of max-separable func- 
tions. Under quite general conditions, explicit solutions for the problem are 
given represented basically in conventional terms rather than in a closed 
form of tropical vector algebra. 

For another problem in IRmax,+ with two-sided equality constraints, 

... t 

minimize p x, 

subject to Ax © b = Cx © d, 

solutions are obtained in [181 ES EH EH]. In particular, [27] presents a 
pseudo-polynomial algorithm to evaluate a solution or to indicate that no 
solutions exist. The algorithm uses an alternating method offered by |37j 
that replaces the equality constraint with two opposite inequalities, which 
are solved by turns iteratively to provide progressively better bounds for a 
solution. A heuristic approach is suggested in [29] to get an approximate 
solution by use of local search techniques combined with iterative procedures 
for solving the problem in particular low-dimensional cases with one and two 
variables. 

2.2 Nonlinear objective functions 

We start with a particular problem that can be written in the form 

minimize (Ax)~d, 
subject to Ax < d. 

The problem is formulated as to find x to provide a best underestimating 
approximation for d by Ax in the Chebyshev norm. It is examined in |15j . 
where a complete explicit solution is given as an application of an abstract 
theory of linear operators over R max ,+ developed there. A similar solution 
is also suggested by [7J. 

Another two problems are initially represented in a different manner, but 
can be put in a form with tropical nonlinear objective function and linear 
constraints. A problem of minimizing a Chebyshev- like distance function un- 
der some constraints is defined in R max ,+ and solved with a polynomial-time 
threshold- type algorithm in [T3]. In fact, with the same algebraic manipu- 
lations as in [23} I3T] , it can be rearranged to take a form with a nonlinear 
objective function and boundary constraints, 

minimize (Ax)~p © q~ Ax, 
subject to g < x < h. 
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A problem of minimization of the range norm is considered in [28] , where 
an explicit solution is suggested in a combined analytic framework that in- 
volves both the semifield K m ax,+ and its dual semifield K m in,+ • The problem 
can actually be written only in terms of lR m ax,+ as 

minimize 

where 1 = (1, . . . , 1) T is a vector of ones (in terms of the semifield R max ,+ )• 
On the basis of a connection between solutions of two-sided equality 
constraints and a mean payoff game established by [38], an iterative com- 
putational scheme is developed in [3§] to find a solution in K max ,+ for a 
problem with nonlinear objective function given by 

■ ■ ■ T i T \— 1 

minimize p x(q x) , 
subject to Ax ffi b < Cx ffi d. 

Now we consider problems that are formulated in terms of a general 
linear ordered radicable semifield and admit explicit solutions in a vector 
form. The analytic technique implemented to solve the problems is based 
on new results in tropical spectral theory and solutions of linear tropical 
inequality developed in [22 | [25 | [1~3]. 

First, a complete solution is given in [40 1 to an unconstrained problem 

minimize x~ Ax © x~p ffi q~x. 

Furthermore, |31] offers a direct explicit solution to another uncon- 
strained problem in the form 

minimize (Ax)~p ffi q~ Ax. 

It is also shown that two problems with linear inequality and equality 
constraints, 

minimize x p ffi q x, minimize x~p(Bq x, 

subject to Ax < x; subject to Ax = x; 

reduce to the former unconstrained problem and so get explicit solutions. 

Finally, in a brief conference paper [20] . a complete direct solution is 
given to a problem 

minimize x~Ax, 
subject to Bx < x 

under the condition that at least one of the matrices A and B is irreducible. 

In this paper, we examine in detail a more general version of the last 
problem with arbitrary matrices and with additional constraints included. 
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3 Preliminary definitions and results 



The purpose of this section is to offer a brief introduction to tropical (idem- 
potent) mathematics so as to provide an appropriate theoretical framework 
for subsequent analysis and solution of tropical optimization problems. 

Both concise introductions to and comprehensive presentations of the 
theory and methods are suggested in a range of published works, including 
[111 151 S3 El El El El [M [ID S31 H21 HU , which sometimes vary in the notation 
system adopted and in the form of results presented. 

Below we give an overview of definitions, notation, and background re- 
sults that mainly follow [25 1 l22| [T3] and form a basis for complete solution 
of the problem under study in an explicit vector form as well as for clear ge- 
ometric illustrations. Additional details and an extensive bibliography can 
be found in the publications listed above. 

3.1 Idempotent semifield 

Let X be a set that is equipped with two operations: addition © and mul- 
tiplication (g>, and contains their related neutral elements: null and iden- 
tity 1. Both operations are assumed to be associative and commutative, 
and multiplication to be distributive over addition. Furthermore, addition 
is idempotent providing x®x = x for all x £ X. Multiplication is invertible, 
which means that there exist a multiplicative inverse x _1 for any x € X + , 
where X + = X\{©} . Endowed with these properties, the algebraic structure 
(X, 0, 1, ffi, (g>) is commonly referred to in the literature as the idempotent 
semifield. 

Due to the associativity of multiplication, the power notation with in- 
teger exponents is introduced in the usual way. For any x € X + and an 
integer p > 1, it is defined as x° = 1, (D p = 0, x p = x p_1 © x = x <g> x p_1 , 
and x~ p = {x~ l ) p . Moreover, we assume that the integer power extends to 
rational exponents, which makes the semifield radicable. 

From here on, we leave out the multiplication sign © for brevity sake, 
whereas the exponents are considered in the sense of the above notation. 

There is a partial order induced on X by idempotent addition so that 
x < y if and only if x © y = y . The definition implies an extremal property 
of addition in the form of inequalities x < x © y and y < x © y . Moreover, 
according to the order, both addition and multiplication prove to be isotone. 

Furthermore, the partial order is assumed to extend to a total order, 
which makes the semifield linearly ordered. Throughout the paper, we in- 
terpret the relation symbols and the problem formulations in terms of this 
linear order. 

Instances of the linearly ordered radicable idempotent semifield include 
Kmax,+ = (K U{-oo}, -oo,0, max, +), R m i n ,+ = (RU {+oo}, +oo, 0, min, +) , 
Kmax.x = (R+U{0},0, l,max, x), and [R miniX = (K+U{+oo}, +oo, 1, min, x), 
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where R is the set of real numbers, R+ = {x G R\x > 0}. 

For example, the semifield R m ax,+ contains the null = — oo and the 
identity 1 = 0. There is an inverse x~ l defined for each x G R and equal 
to —x in conventional algebra. For any x,y G R, the power x y coincides 
with the arithmetic product xy . The order induced by idempotent addition 
corresponds to the natural linear order on IR . 

The semifield K m in,x is equipped with (D = +oo and 1 = 1. The in- 
verse and power notations have the standard interpretation. The relation 
< defines a reverse order to the linear order on IR . 

Finally note that all above semifields are isomorphic to each other. 

3.2 Idempotent semimodule 

Consider X n to be a Cartesian power of X with elements taken to be column 
vectors. There is a vector (D G X n that has all components equal to . Any 
vector without zero components is called regular. The set of all regular 
vectors in X n is represented by X™ . 

For any vectors a = (aj) and b = (6j) from X ra , and a scalar x G 
X, vector addition and scalar multiplication are defined component-wise 
according to the usual rules 

{a © b}i = en © {xa}i = xai. 

Endowed with these operations, the Cartesian power X n forms an idem- 
potent semimodule over X . 

The extremal property of idempotent addition extends in the semimod- 
ule to component-wise vector inequalities a < a © 6 and b < a © b. The 
isotone properties of scalar operations imply that both vector addition and 
scalar multiplication are also isotone in each argument. 

To provide a clear visualization of key phenomena inherent in idempo- 
tent semimodules, we use graphical illustrations in terms of the semimodule 
Jsmax , . Any vector in IR^ax + is shown on the plane with a Cartesian coor- 
dinate system that can be supplemented by two artificial points if needed. 
The points are respectively put before the left end of the horizontal axis and 
below the bottom of the vertical axis to indicate, in a schematic way, the 
value = — oo required to represent zero coordinates of irregular vectors. 

A geometric interpretation of the operations in the context of the semi- 
module IRmax + is given in Figures [T] and [2j Idempotent addition of two 
vectors with a common base point at the origin of the coordinate system 
goes by a "rectangle rule" . According to the rule the sum is obtained as the 
upper right vertex of an upright finite (Figure [TJ left) or infinite (Figure HJ 
right) rectangle that has vertices at the end points of the vectors to add. 

Scalar multiplication of a vector by a nonzero scalar results in moving 
the end point of the vector along a line at 45° to the coordinate axes for 
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a © 6 



61 ai 




b\ = —00 ai 



Figure 1: Vector addition of regular vectors (left) and of regular and irregular 
vectors (right). 




xa 
V 



a 



xa\ 




Figure 2: Scalar multiplication of a regular vector (left) and of an irregular 
vector (right). 



regular vectors (Figure [21 left) or along an artificial vertical (horizontal) line 
for irregular vectors (Figure [21 right). 

As usual, a vector b is said to be linearly dependent on vectors 01, . . . , a m 
if b = x\a\ © • • • © x m a m for some scalers x\ , . . . , x m . Specifically, two vec- 
tors a and fe are collinear if there is a scalar x such that fo = xa. 

Given vectors 01, . . . , a m , the linear combinations x\a\ © ■ ■ ■ ffix m a m for 
all x\, . . . , x m form a linear span of the vectors and constitute an idempotent 
subsemimodule. Graphical examples of the linear span in R^ax + are snown 
in Figure as areas surrounded by hatched borders. 

The linear span of two regular vectors is a strip bounded by the lines 
drawn through the end points of the vectors (Figure O left). When one 
of the vectors is irregular, the strip extends to take a form of a half-plane 
(Figure El right). 

For any nonzero column vector a = (ai) , we define its multiplicative 
conjugate transpose as a row vector a = (a^) with components aj = a^ 1 
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Figure 3: Linear span of regular vectors (left) and of regular and irregular 
vectors (right). 

if a,i > , and a~ = otherwise. If both vectors a and b are regular, then 
the component- wise inequality a < b implies a>b. 

3.3 Matrix algebra 

Consider matrices with entries from X. For conforming matrices A = (ay), 
B = (hj), and C = (c^), and scalar x € X, matrix addition and multipli- 
cation together with scalar multiplication follow the element-wise formulas 

{A © B}ij = dij © b^ , {AC}ij = OjfcCfcj, {xAjjj = xa^. 

k 

Specifically, multiplication of a matrix A = (a^) by a vector a: = (x«) 
of order n gives a vector with components 

= (ii\Xi © • • • © di n x n . 

The matrix operations are element- wise isotone in each argument. For 
any matrices A and B of the same size, it holds A < A © B and B < 
A®B. 

A matrix that has only zero entries is the zero matrix represented by . 

Consider square matrices of order n over X and denote the set of the 
matrices by X raxra . Any matrix having all off-diagonal entries equal to 
is a diagonal matrix. A diagonal matrix with 1 along the diagonal is the 
identity matrix denoted by I . 

The matrix power is introduced in a regular way. For any matrix A € 
X nxn and integer p > 1, it is defined as A = I , A p = A p ~ l A = AA p ~ l . 

If all entries of a matrix above or below the diagonal are zero, the ma- 
trix is triangular. A triangular matrix with zero diagonal entries is strictly 
triangular. 
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A matrix is reducible if simultaneous permutations of rows and columns 
can put it into a block-triangular normal form, and irreducible otherwise. 
The lower triangular normal form of a matrix A is given by 



/A n ... \ 
A 2 i A22 



(1) 



\ A s i A s2 . . . A ss ) 

where An is either irreducible or zero matrix of order m , whereas A^ 
is an arbitrary matrix of size rij x rij for all i = l,...,s, j < i, and 

ni H h n s = n. 

The trace of any matrix A = (o«) is routinely defined as 

tr A = an © • • • © «nn- 

For any matrices A and £?, and scalar x it holds 

tr(A©S) =trA0trS, tr(AS) = tv(BA), tr(xA) = xtr(A). 

Application of these properties leads, in particular, to a binomial identity 
for traces derived by |20| in the form 

m 

tr(A0 J B) m = tr ( AS* 1 • • • AB ik ) © tr B m , (2) 

fc=l iiH hifc=m— fc 

which is valid for any non-negative integer m. 

As usual, a scalar A is an eigenvalue of a matrix A, if there exists a 
nonzero vector x such that 

Ax = Xx. 

Every irreducible matrix has only one eigenvalue, whereas reducible ma- 
trices may possess several eigenvalues. The maximal eigenvalue (in the sense 
of the order on X ) is called the spectral radius of A and directly calculated 
as 

n 

X = 0tr 1 / m (A m ). 

m=l 

The spectral radius A of any matrix A offers a useful extremal property 
considered by [2T], which holds that 

min x Ax = A, 

where the minimum is taken over all regular vectors x. 
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4 Explicit solution to a linear inequality 



Explicit results for optimization problems to be derived in the next section 
make well use of direct solutions to the following problem. Given a matrix 
A € X nx ™ and vector b G X n , find all vectors x £ X™ to satisfy the 
inequality 

Ax®b<x. (3) 

There is an equation with the same terms as in ([3j), which is frequently 
referred to as the non-homogeneous Bellman equation. The equation is 
examined in a number of works, including [HI [6l SSI [12j HQ El |9] , whereas 
solution to the inequality seems to receive less attention in the literature. 

Complete solutions for the inequality with both irreducible and reducible 
matrices can be found in |25l IT3] , where they are obtained by reducing the 
inequality to an equation. Below we concentrate on regular solutions in the 
case of general reducible matrices and offer a direct proof to get an explicit 
solution in a slightly different and more compact form. 

For any matrix A, we define a function that is given by 

Tr(A) = trA0---©trA n , 
and a star operator that takes A to a matrix 

a* = /0 Ae---eA n_1 . 

Both the function and the star operator act together in an inequality for 
matrix powers that represents, in somewhat different terms, a known result 
apparently first suggested by [44] . The statement, which is referred to below 
as the Carre inequality, holds that the condition Tr(A) < 1 involves, for all 
integer k > 0, the inequality 

A k < A*. 

For irreducible matrices, a solution to ([3]) can be represented as follows. 

Lemma 1. Let x be the general regular solution of inequality Q with an 
irreducible matrix A . Then the following statements hold: 

1. If Tr(A) < 1, then x = A*u for all regular vectors u >b. 

2. If Tr(A) > 1, then there is no regular solution. 

Figure H] offers examples of solutions to inequalities with a common ma- 
trix A = (04,0.2) and different vectors b in IR^ ax + . The set of solutions is 
indicated by semi-infinite regions with hatched borders. Together with the 
columns of the matrix A, those of A* = (a^Og) are also presented. 

In the rest of the section, we show that the result of LemmaUJcan directly 
be extended to solve the inequality with irreducible matrices. We start with 
an auxiliary result that provides a useful representation for (A © B)* and 
can be considered as an analogue of an identity that appears in [6j. 
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Figure 4: Solution to linear inequalities with an irreducible matrix. 

Proposition 1. For any matrices A and B with Tr( A © B) <t, it holds 

(A © B)* = (A*B)*A*. (4) 

Proof. To verify the identity, it is sufficient to show that two opposite in- 
equalities between its sides are valid. By expanding both left and right parts 
of (JU) and rearranging the terms, we obtain one of the inequalities 

n-1 

(AffiS)* = A i °(BA il ) ■ ■ ■ (BA im ) 

m=0 0<io+iiH \-i m <n—m—l 

n-1 

^0 A ia (BA il ) ■ ■ ■ (BA im ) = (A*B)*A*. 

m=0 0<«o,!i,...,im<n- 1 

Furthermore, we denote C = A © B and see that A < C and B < C . 
Application of the Carre inequality gives 

n 2 -l 

{A*B)*A* < (C*C)*C* = C m = C* = {A © B)*, 

m=0 

which shows that the opposite inequality is also true. □ 

Now we are in a position to present the main result of the section. 

Theorem 2. Let x be the general regular solution of inequality ([3]) with an 
arbitrary matrix A . Then the following statements hold: 

1. If Tr(A) < 1, then x = A*u for all regular vectors u >b. 

2. If Tr(A) > t, then there is no regular solution. 
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Proof. If the matrix A is irreducible, then the statement of the theorem is 
provided by Lemma [TJ Suppose that A is reducible. Provided that A = 0, 
the theorem is trivially true. Consider the case when A / C and assume, 
without loss of generality, the matrix A to have the form of JT]). 

We define block-diagonal and strictly lower block-triangular matrices 



D 



( An 




1'2 1 



\A sl ... A 



s.s—l 



\ 



/ 



to be corresponding parts in the additive decomposition A = D © T. 

According to the row partition of the matrix A in the form (fTJ) , we write 
both vectors x and b in the block form 



( X1 \ 


b = 


( bl \ 


V x s J 







X 



where Xi and bi are vectors of order for each % = 1, . . . , s. 

Furthermore, we represent ([2D as a system of inequalities, one for each 
element of the partition, 

A{\X\ © • • • © AaXi © bi < Xi, i = 1, . . . , s. 

Suppose that Tr(A) < 1. In this case, for each i, we have 

Tr(A rt ) < TV(Au) © • • • © Tv(A ss ) = Tr(A) < 1. 

It follows from Lemma [T] that inequality i can be solved with respect to 
Xi . Assuming Vi to be a regular vector of order , the solution is given by 



Xi = A* { Vi, Vi > Aaxt 



Aj j_ \Xi © bi. 



With another notation Ui to represent a regular vector of order rii , the 
last inequality can be rewritten as 

Vi = A a xi © • • • © A iyi -iXi © Ui, Ui > hi. 

Turning back to the solution of inequality i , we arrive at the represen- 
tation 



Xi — AaiyAi\X\ 



© Ai i-iXi) © A*uUi, Ui>bi. 



After concatenation of all vectors u\ , . . . , u s into one vector u of order 
n , we write all solutions in the form of an implicit equation for x , 



x = D*Tx © D*u, 
14 



u> b. 



a 2 




a 2 




Figure 5: Solution to linear inequalities with a reducible matrix. 

Considering a strictly block-triangular form of T and a conforming 
block-diagonal form of D* , it is not difficult to verify that (D*T) r = 
for all r > s. 

By iterating the implicit equation, we arrive at a regular solution 



x 



(D*T) 



© (D*T) a ~ 1 )D*u = (D*T)*D*u, u>b. 



With application of identity to A = D © T , we put the solution in 
the final form 

x = A*u, u > b. 

If Tr(A) > 1 then there is at least one row partition i with Tr(An) > 1. 
In this case, according to Lemma [H no regular solutions X\ exist for the 
partition and so no regular solutions x for ([3]). □ 

A graphical illustration of solutions to the inequality with a common 
reducible matrix A in Rj^ax + ^ s gi ven i n Figure [5j 



5 A constrained optimization problem 

We start with a real-world problem drawn from project scheduling and then 
represent it in terms of the semifield R m ax,+ • The problem serves as mo- 
tivation and illustration for the solution of a more common problem that 
is formulated in the context of a general idempotent semifield, which cov- 
ers [R m ax,+ as a specific instance. As the main result, we obtain a solution 
to the common problem under fairly general assumptions. Particular cases 
of the problem are considered and illustrated with numerical and graphical 
examples in the framework of the semifield K m ax,+ • 
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5.1 A project scheduling problem 

Below we concern with optimal scheduling of collection of activities (jobs, 
tasks), which presents a rather typical problem of schedule development in 
project management (see, e.g., [47, 48J for further details and examples). 

Consider a project that involves n activities operating under precedence 
constraints, which are referred to as Start-to-Finish, Start-to-Start, and 
Early Start constraints. The Start-to-Finish constraints do not allow an 
activity to complete until predefined times after initiation of other activi- 
ties. The activities are assumed to complete as early as possible to meet 
the constraints. The Start-to-Start constraints determine minimum allowed 
times between initiation of activities, whereas the Early Start constraints 
specify earliest possible dates for activity initiation. 

Each activity in the project is characterized by its flow time (also known 
as turnaround or processing time) defined as the time interval between its ini- 
tiation and completion. The problem is to design a schedule that minimizes 
the maximum flow time over all activities, subject to the above constraints. 

For each activity i = 1, ... ,n, we denote its initiation time by X{ and 
its completion time by yi . Let aij be a minimum possible time lag between 
initiation of activity j = 1, . . . ,n and completion of i. 

Given , the completion time of activity i , due to the Start-to-Finish 
constraints, must satisfy the relations 

Xj + aij<yi, j = l,...,n, 

where at least one inequality holds as an equality. Note that, in the case 
when ctij is not given for some j , it is assumed to be — oo. 
The relations can be combined into one equality of the form 

yi = max(xi + an,...,x n + a in ). 

By replacing the ordinary operations with those in IR max ,+ , we get 

Vi = Q-ji^i ® • • • © cii n x n , i = 1, . . . , n. 



With a matrix and vectors, 





( an . 


■ ain \ 




( *1 ^ 




( m \ 


A = 






X = 




1. 








ann / 








\Vn ) 



we represent the scalar equalities in a vector form 

y = Ax. 
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Furthermore, we formulate a problem of minimizing the maximum flow 
time over all activities in the project. In ordinary notation, the objective 
function in the problem is given by 

max(yi -x 1 ,...,y n - x n ). 

Since, in terms of R maXj+ , the function is written as 

x^yi ■ ■ ■ ® x~ l y n = x y = x Ax, 

we arrive at a problem to find all regular x so as to 

minimize x Ax. 

Now we incorporate Start-to-Start and Early Start constraints into the 
problem formulation. For each activity i = l,...,n, we denote by gi an 
early possible initiation time and by Cjj a minimum possible time lag be- 
tween initiation of activity j = 1, . . . ,n and initiation of i. Given Cjj and 
gi , the initiation time Xi for activity i is subject to the relations 

max(xi + Cji , . . . , x n + Q n ) < Xi, 

9i < Xi. 

Representation in terms of R m ax,+ gives scalar inequalities 

9i < Xi. 

With the matrix-vector notation C = (c^) and g = (gi) , we get 

Cx < x, 
g < x. 

The optimal scheduling problem under consideration takes the form 

minimize x Ax, 
subject to Cx < x, 
g < x. 

We note in conclusion that the solution of the problem (if exists) is not 
unique. Specifically, if there is a solution x , then ax , where a > , is also 
a solution. In the context of project scheduling, however, a unique solution 
naturally appears when the minimal solution is of interest. 



17 



5.2 A general representation 



Let X be a linearly ordered radicable idempotent semifield, A,C £ x nxn 
be given matrices and g € X n be a given vector. The problem is to find all 
regular solutions x £ X™ so as to 

minimize x~Ax, 

subject to Cx < x, (5) 
g <x. 

Since, for all vectors x, the inequalities Cx < x and x > g can readily 
be combined into one equivalent inequality Cx © g < x, the problem has 
an alternative representation in the form 

minimize x Ax, 

(6) 

subject to Cx © g < x. 

Consider the set of regular vectors that are determined by the inequality 
constraints in the problem and find out when the set is nonempty. Clearly, 
an appropriate necessary and sufficient condition follows from Theorem [2] 
and requires that Tr(C) < 1. 

5.3 The main result 

A complete solution to the optimization problem is given as follows. 

Theorem 3. Suppose A is a matrix with spectral radius A > , C is a 
matrix with Tr(C) < 1, and g is a vector. Define a scalar 

n-1 

= A + tr 1 ^ (AC* 1 •••AO) (7) 

k=l 1<«iH yi^<n—k 

and a matrix 

S = d- 1 A®C. (8) 
Then the minimum in ^ is equal to 6 and attained if and only if 

x = S*u, u>g. (9) 

Proof. We start with the problem in the form of jB]) and then follow the 
approach suggested in [20]. We introduce an additional variable so as to 
reduce the problem to an inequality. Furthermore, existence conditions for 
solutions of the inequality is used to determine the variable, whereas the 
solution of the inequality is taken as the solution of the problem. 
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Suppose 9 is the minimum of the objective function in the problem and 
note that 9 > A > 0. All regular vectors x that yield the minimum are 
therefore determined by the system 

x~Ax = 9, 
Cx®g <x. 

Consider the first equality x~ Ax = 9. Since for all x we actually 
have x~ Ax > 9, the equality can be further replaced with an inequality 
x Ax < 9. 

Next we take the inequality x~ Ax < 9 and multiply it by 9 _1 x from the 
left. Since xx > I for all regular x, we have 9~ l Ax < 9~ 1 xx~ Ax < x 
and so get a new inequality 9~ l Ax < x. 

Considering that left multiplication of the last inequality by 9x~ and 
an obvious equality x~x = 1 give the former one, both inequalities are 
equivalent. 

Now the solution set of the problem is given by the system of inequalities 

6~ x Ax < x, 
Cx®g <x. 

It is easy to verify that the system is equivalent to one inequality, 

Sx®g<x, (10) 

where the matrix S is defined as in ([8]). 

It follows from Theorem [2] that a necessary and sufficient condition for 
regular solutions of the inequality to exist is given by 

Tr(S) = Tr^A ffl C) < 1. (11) 

To further examine the existence condition, we rewrite the left side of 
the inequality and then apply binomial identity ([2]) to get 

rt 

Tr(5) = Qtrr^eC)™ 

m=l 

n m 

= 00 9~ k tr(AC il ■ ■ ■ AC ik ) © Tr(C). 

m=l k=l iiH Mfc=m— k 

Furthermore, rearrangement of terms in the last sum gives the following 
expression 

n 

Tr(5) = 9- k tr ( AC h ■ ■ ■ AC ik ) ©Tr(C). 

k=l 0<iiH hifc<n— k 
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Considering that Tr(C) < 1 by the conditions of the theorem, existence 
condition ([TT]) reduces to inequalities 

0~ k tv(AC il ■ ■ ■ AC ik ) < 1, k = l,...,n. 

0<iiH h*fc<«— k 

By solving the inequalities with respect to 9, we arrive at inequalities 
9> ti 1/k {AC h ■■■ AC ik ), k = l,...,n, 

0<iiH \-ik<n~ k 

or, equivalently, at one inequality 

n 

^-0 tr l / k {AC n ••• AC ik ) 

k=l 0<iiH \-ik<n—k 

= X ® tr 1/k ( ACh ■■■AC* k ). 

k=l l<iiH \-i^<n—k 

In order for 9 to be the minimum in the problem, the last inequality 
must be satisfied as an equality, which gives 0. 

Finally, application of Theorem [2] to inequality (TlOj) leads to the solution 
in the form of ([9]). □ 

Note that the minimum solution of the problem is given by xq = S*g . 
5.4 Particular cases 

First assume that C = and consider a problem 

minimize x~Ax, 

(12) 

subject to x > g. 

The next solution to the problem is a direct consequence of Theorem [3l 

Corollary 4. Suppose A is a matrix with spectral radius A > and g is a 
vector. Then the minimum in (|12p is equal to A and attained if and only if 

x = (\~ 1 A)*u, u>g. 

Now consider the case when a problem has no lower bound constraints 
and so is given by 

minimize x Ax, 

(13) 

subject to Cx < x. 
By setting g = in Theorem [3j we immediately get the following result. 
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Corollary 5. Suppose A is a matrix with spectral radius A > and C is 
a matrix with Tr(C) < 1. Define a scalar 6 by ^ and a matrix S by ([8]). 
Then the minimum in (|13p is equal to 6 and attained if and only if 



x = S*u, u > 0. 

It is easy to see that, when at least one of the matrices A and C is 
irreducible, the last result coincides with that in |20j. 

Finally, for a problem without any linear constraint in the form 

minimize x~Ax, (14) 

the solution offered by the theorem reduces to the following. 

Corollary 6. For any matrix A with spectral radius A > 0, the minimum 
in (fT4"j) is equal to A and attained if and only if 



x = (A 1 A)*u, u > 0. 
Note that this result is quite consistent with that in [40J. 

5.5 Illustrative examples 

To illustrate the results obtained for problem ([5]) with both irreducible and 

v max,+ • 



reducible matrices, we give numerical and graphical examples in R 2 



5.5.1 Irreducible matrices 

We start with a problem with irreducible matrices A and C, which are 
defined, together with a vector g, as follows: 

-2 A / -10 \ _ / -9 

-7 -3 J ' ° " I 4 -3 r 9 ~\ 6 



First consider problem (|14[) . which has no constraints. Taking into ac- 
count the matrices A and 

4 o / -2 \ / (.) -2 \ ( -2 
A 



-7 -3 y V -7 -3 J \ -7 -6 
we find the spectral radius A = tr A © tr 1//2 (A 2 ) = = 1. Then we get 

(A-U)* = A* = I(BA= ( _° 7 ~ 2 Q 

The solution to the unconstrained problem takes the form 

Q ^ 

x - I _ 7 ) «, we c " 
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which is graphically represented as in Figured] (left). 

Furthermore, we obtain the solution to the inequality constraints in the 
form Cx © g < x. To apply Theorem [21 we take the matrices C and 

C 2 = 

\ 4 -J y \ 4 -6 , 
Then we verify that Tr(C) = tr C © tr C 2 = = 1. After calculating 

-10 








J 




" J 


-10 \ 


4 


--3 




--3 




-6 J 



C: - J ■ C ■ ! () 

we arrive at the solution (see Figure EJ right) 

/ -10 \ / -9 

X= [i J ^( 6 

Now we use Theorem [3] to handle the constrained problem. We calculate 

. „ / -2 \ / -10 \ / 2 -5 
AC ={-7 -3 JU -3 ) = {l -6 

and then get 9 = A © tr( AC) = 2. Furthermore, we have 

s = e-i A ec=(l s* = ies=( A ~ 4 Q 

All solutions to the problem take the form 

x= { \ "o) u ' U ~{~1 

To simplify the solution, note that the columns in the matrix S* are 
collinear. Therefore, we can write 

x - $' u - ( ° 4 ~ Q J u = ( ° ) ( - I ) w. 



With a new scalar variable 

v = ( -4 ) « > ( -4 ) flf = ( -4 
the solution of problem ([5]) gets its final form 

x= ( 4 ) U ' U - 2 ' 

Figure (right) combines the solutions constraints with the final solution 
of the entire constrained problem. The final solution is depicted with a thick 
half-line that coincides with the right border of the feasible area defined 
by the constraints. The half-line starts at the end point of the vector xq 
representing the minimal solution. 
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Figure 6: Solutions to unconstrained (left) and constrained (right) problems 
with irreducible matrices. 

5.5.2 Reducible matrices 

Now we solve problem ([5]) with reducible matrices under the conditions that 

"=(-»=!)• ■=(!)■ 

We successively calculate 

* = (zt , a = o, ( A-u r = A. = ( J --). 

The solution to the problem without constraints is given by (see Figure El 
left) 

f -oo\ n2 
x = { _ 4 J " € R ■ 

Furthermore, we obtain the solution to the inequality constraints. Since 

°"=(-»:S)- ^°>=* c -(-oo°"o). 

the solution takes the form (FugureO right) 

~U°D- ■*(:)■ 

To solve the constrained problem, we apply Theorem We get 
and then write the solution of problem ([5]) in the form 
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Figure 7: Solutions to unconstrained (left) and constrained (right) problems 
with reducible matrices. 

The solution is shown in Figure[7] (right) as a half-strip region surrounded 
by a thick hatched border. 

6 Conclusions 

We started this paper with an overview of known tropical optimization prob- 
lems, which demonstrates that tropical optimization is currently under in- 
tense development as a promising solution approach for many real-world 
problems in operations research. The paper examines a new optimization 
problem with a nonlinear objective function and linear inequality constraints 
in a rather general setting. The problem sufficiently extends that in [20] by 
eliminating irreducibility conditions for matrices involved as well as by in- 
troducing additional constraints. We offer an example problem drawn from 
project scheduling as a real motivation and clear illustration of the study. 

The derivation of the solution is based on definitions, notations, and basic 
results developed in \25\ [22j [T3] , which are outlined in a brief introduction 
supplemented by appropriate graphical illustrations. As a preliminary step, 
a new complete solution to a linear inequality was obtained for the case of 
arbitrary matrices. The solution extends known results in [25} [T3] and is of 
independent interest. 

The solution to the optimization problem under study follows the ap- 
proach suggested in [20] . which reduces the problem to the solution of a 
linear inequality with a parametrized matrix. By applying the above result 
for linear inequalities, we obtain a complete solution to the optimization 
problem in a compact vector form. To illustrate the result, numerical solu- 
tions and graphical representations are given for two-dimensional problems. 

Future research is expected to concentrate on further extension of the 
problem intended to account for new types of constraints, including equality 
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constraints, and more general form of the objective function. Another line 
of investigation will be development of new real-world applications of the 
results. 
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